Longitudinal Fluctuations in Partonic and Hadronic Initial State 
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Collective flow in collisions between Lead nuclei at LHC are influenced by random initial state 
fluctuations, especially for odd harmonics. Here we extend fluctuation studies to longitudinal fluctu- 
ations, which may have significant effect on the rapidity distribution of odd harmonics. Furthermore 
center of mass rapidity fluctuations are measurable, but not yet analysed. Here in the PACIAE par- 
ton and hadron molecular dynamics model we make an analysis of initial state fluctuations. As 
previous analyses discussed mainly the effects of fluctuations on eccentricity and the elliptic flow we 
pay particular attention to the fluctuations of the Center of Mass rapidity of the system, which is 
conservatively estimated in our model as AycM = 0.1, by neglecting all pre-equilibrium emission 
effects that are increasing the ycM fluctuations. 

PACS numbers: 12.38.Mh, 25.75.-q, 25.75.Nq, 51.20.+d 



I. INTRODUCTION 

Global collective observables are becoming the most 
essential observables in ultra-relativistic heavy ion 
reactions [1]. When we want to extract precise knowl- 
edge from experiments, both on the Equation of State 
(EoS) and the transport properties of matter 0, [1] we 
have to invoke a most realistic description with fully 3+1 
dimensional dynamical evolution at all stages of the reac- 
tion, including the initial state. This most adequate de- 
scription of all stages can only be achieved by the multi- 
module, hybrid models. (See e.g. refs. Q.) 

The initial state, where we have very little direct ex- 
perimental information, is of paramount importance in 
the theoretical description. This leads to a wide variety 
of initial state models, which behave differently. Theo- 
retical models and experimental results indicate that the 
initial state fluctuations are essential in understanding 
the data, although in the global continuum (fluid dynam- 
ical or field theoretical) models these fluctuation effects 
may inherently not be present and even may not survive 
to the hadronic final state. Nevertheless, we need to an- 
alyze the behavior of these initial state models from the 
point of view of fluctuations. (See e.g. refs. Q.) 

However, one has to take into account that the Cen- 
ter of Mass (CM) rapidity is not exactly the same for all 
events because of random fluctuations in the initial state 
caused by the difference of participant nucleon numbers 
from projectile and target. This leads to considerable 
fluctuations at large impact parameters, where the flow 
asymmetry is the strongest, but the number of partici- 
pant nucleons is the smallest. 

Just as all initial state fluctuations, we have two 
sources of CM-rapidity fluctuations: First, the number 
of nucleons are randomly located in the configuration 
space and due to their fluctuating location, the num- 
ber of participants from the target and projectile nucleus 
must not be the same event- by- event, even in the sym- 



metric, A+A, collisions. Second, those nucleons, which 
are in the geometrical participant zone, may actually not 
collide with any single nucleon from the opposite nucleus, 
consequently these will not become participants. Some 
recent results on the subject concerning the V2 and i>3 
fluctuations are discussed in refs. 



Up to now less attention is paid to the fluctuations 
in the beam direction. The expected momentum and/or 
rapidity fluctuations in this direction may be bigger due 
to the large beam momentum in recent experiments. In 
case of CM-rapidity fluctuations there is an additional 
problem: It is not obvious how tightly bound system 
is the initial state. The number of participant nucleons 
may not come from the projectile and the target nuclei 
equally, there can be one or a few more nucleons from 
one side. The momentum carried by the extra nucleons, 
may be shared (i) by all participants equally in a tightly 
bound system (a single large confined QGP bag, may be 
considered as such a system) or (ii) by a loosely connected 
cloud of nucleons (where the extra nucleons have little di- 
rect effect on the participant matter). In the later case, 
although the total momentum is conserved, the internal 
energy of the participant matter is increased considerably 
by the energy of the extra nucleons but the momentum 
of the participant matter is not correlated with the mo- 
menta of the extra nucleons. So, the collective rapidity 
change is much less. 

It is important to mention that the phase transitions 
and the consequent fluctuations both in and out of QGP 
may enhance the collective behavior of the system [9[. 
However, it is rather difficult to estimate the conse- 
quences of such transitions and fluctuations to the CM- 
rapidity fluctuations. From the point of view of initial 
state fluctuations we have to arrive at system, which is 
close to local equilibrium, thus, at high energies the tran- 
sition to QGP has to happen earlier than the formation 
of the initial state. Thus, it is important to study the 
CM-rapidity fluctuation as an observable on its own to 
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learn about energy deposition, and also due to its strong 
effect on flow observables. (See e.g. ref. [13]). 

In this work, after some simple considerations, we 
present an analysis of these fluctuations in the PACIAE 
model, where the major sources of fluctuations are taken 
into account. 



II. ANALYTICAL ESTIMATES FOR THE 
CM-RAPIDITY FLUCTUATIONS 



As mentioned above, the initial state fluctuation is 
stemming from the participant nucleon number (Na + 
Nb = N part ) fluctuation. Here Na and Nb are the num- 
bers of participant nucleons from the projectile and tar- 
get nuclei, respectively. The participant matter forms 
then the initial state system. In the following examples 
we present three situations where different fraction of the 
beam energy is contributing to the total transverse mass 
of the locally equilibrated participant matter. 

Let us first estimate the effect of fluctuations of the 
participant matter for a impact parameter of b = 0.7b max 
collision in Pb+Pb reactions at the LHC energy of 
1.38 + 1.38 A-TeV, for a tightly bound and unexcited 
system. We assume that one extra nucleon from the pro- 
jectile nucleus will be absorbed into the participant mat- 
ter, which otherwise would contain N part = Na+Nb = 
32.5+32.5 = 65 nucleons. Then this extra projectile nu- 
cleon, 5N = Na—Nb = 1, carries m t * sinh(yo) momen- 
tum, where yo = 8 is the beam rapidity at the above 
LHC energy and mt = ton is the transverse mass of a 
nucleon in the beam. If this extra momentum is absorbed 
in the participant matter, then according to the momen- 
tum conservations: 

P z = M t CM smh(A 2 , CM ) = SNm t smh(y ) (1) 
E = M t CM cosh( Aj/cm) = N part m t cosh(yo) (2) 

this extra nucleon will lead to a change of the CM- 
rapidity, AycM (which is zero if the participant nucleons 
are coming in equal numbers from the projectile and tar- 
get) . In the above equations the M^ M is transverse mass 
of the participant matter. 

In the initial state model based on expanding flux 
tubes or streaks ll| used in fluid dynamical calculations 
[Tol HI], the initial state system is tightly bound and 
stopped within each "streak" . Thus, this model is appli- 
cable streak by streak and its momentum change is more 
pronounced for the peripheral streaks where the asym- 
metry between the projectile and target involvements is 
the biggest. In this initial state model the transverse 
mass, 

M CM 

is more than what would arise from the 
nucleon masses, N par tmN, due to the field strength in 



rapidity can be derived from Eq. (p}: 



the string. So M t ° = N part (m t + La), where L is 
the length of the streak and a is the effective string 
tension. If the participant matter is weakly excited, 
Mf M w N par t {m t + IGeV). The resulting shift of CM- 



AycAf arsinh 



SN m t 



N P art{m t + IGeV) 



sinh(yo) 



= 3.1. 



Thus, CM-rapidity fluctuations may be quite substantial. 
In this case a large fraction of beam energy should be 
carried away through other channels, like pre-equilibrium 
emission. 

For the initial state in hadronic transport models the 
momentum of extra nucleons are hardly influencing the 
momenta of the other participant nucleons. The extra 
nucleons are not stopped in this picture, the transverse 
mass (Mf M ) in the above expression includes large pre- 
thermal momenta, but M t can still be proportional to 
m t * sinh(y ). In such a model the CM-rapidity fluctu- 
ation will be significantly smaller. For example, in the 
above b = 0.7b max Pb+Pb reaction at (1.38+1.38) A-TeV 
if we assume 65 + SN, (where 5N = 1) participant nu- 
cleons and full equilibration, so that 2 /3 rd of the beam 
kinetic energy is converted into the transverse mass of 
the participant matter, and Mf M can be approximated 
as Mf M = N par t(mt + e * 2/3) where e = 1.38TeV per 
nucleon in the Lab/CM frame. Then the CM-rapidity 
fluctuation can be approximated as 



A ycM ~ arsinh 



SN m t 



_Npart(m t +2eo/3) 



sinh(yo) 



0.025. 



(3) 

Although here we discuss the hadronic initial state in a 
hadronic transport model, it is suitable for the partonic 
initial state in hadron and parton transport models also. 

The other limiting case is when all reaction energy is 
absorbed in the participant matter. Then both Eqs. (|1I2[) 
are satisfied, and for the same example of Pb+Pb colli- 
sion as above the resulting CM-rapidity is 



A yen — artanh 



5N 



N 



part 



tanh(yo) 



= 0.015 



(4) 



The above considerations show that the question of 
initial state fluctuations is a rather complex and model 
dependent question. After all, the collectivity or loose- 
ness of the initial state must be estimated experimentally. 
The CM-rapidity fluctuations may provide a very good 
tool to this research. 



III. LONGITUDINAL FLUCTUATIONS IN 
PARTONIC INITIAL STATE IN PACIAE MODEL 

We discussed above the hadronic initial state, now we 
turn to the partonic initial state. In the parton and 
hadron cascade model, PACIAE Jl3| the initial partonic 
state is generated as follows: 

1. The overlap zone and the number of participant nu- 
cleons from the projectile and target are first cal- 
culated geometrically [l4| for an A+A (or A+B) 
collision, at a given impact parameter. 
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2. The participant nucleons are distributed randomly 
inside the overlap zone, starting from nucleons in- 
side the corresponding nuclear sphere having an 
isotropic Woods-Saxon distribution. Nucleons are 
given beam momentum, and a particle-list of initial 
nucleons is constructed. 

3. An A+A (A+B) collision is decomposed into 
nucleon-nucleon (NN) collision pairs and every one 
with a collision time calculated by assuming that 
the nucleons propagate along straight line trajecto- 
ries and interact with the NN inelastic (total) cross 
sections. Then the initial NN collision-list is con- 
structed by these NN collision pairs. 

The PACIAE model assumes that if a NN colli- 
sion happens both colliding nucleons become par- 
ticipants, and eventual occupations of final parti- 
cle states are disregarded. These approximations 
would decrease the longitudinal fluctuations and 
angular asymmetries [15| . 

4. A NN collision pair with the earliest collision time 
is selected from the collision list, and the final state 
of the collision is obtained by the PYTHIA model 
with string fragmentation switched-off. Afterwards 
the diquarks (anti-diquarks) are broken randomly 
into quark pairs (anti-quark pairs), and one obtains 
a configuration of quarks, anti-quarks, and gluons, 
beside a few hadronic remnants for a NN collision. 

Although gluons are treated as point like parti- 
cles, this treatment is not accurate, as gluons are 
mediating the interaction among the color charges 
and they have significant role in the formation and 
hadronization of QGP. In these transitions, the en- 
ergy of gluons is connected to the masses of the 
hadrons and to the energy of the emitted high 
energy photons. We neglect photons in the ini- 
tial state and so we neglect pre-equilibrium photon 
emission also. The detailed treatment of the gluons, 
hadron-parton transition and pre-equilibrium emis- 
sion would increase the ycM fluctuations. To in- 
clude these effects would be overly complicated and 
not realized in models similar to PACIAE. Thus, in- 
stead we chose to neglect the gluon contribution to 
Vcm- In the present highly approximate treatment, 
where gluons are treated as point like classical par- 
ticles, the inclusion of the gluons would reduce ycM 
fluctuations contrary to the physical expectations. 

5. Each of the particles (nucleons) travels along 
straight line trajectories between two consecutive 
NN collisions. After the collision the particle list 
and collision time list are updated, the last step 
and this process are repeated until the NN colli- 
sion list becomes empty (the NN collision pairs are 
exhausted) . 

The hadron and parton cascade model, PACIAE, in- 
cludes the most important geometrical effect of the fluc- 
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FIG. 1: Initial state fluctuation of the number of extra nu- 
cleons, 6N, in 100+100 A-GeV 0-5% central and 70-80% pe- 
ripheral Au+Au collisions in PACIAE model. 



tuation of center of mass momentum in heavy ion colli- 
sions, as the positions of the initial nucleons are random 
following the original Woods-Saxon profiles of the pro- 
jectile and target nuclei. Then in the overlap region nu- 
cleons may collide with each other according to NN cross 
section, and those which do not, will become spectators. 
This construction provides the participant nucleons, their 
positions and momenta, as well as the number of spec- 
tators from the projectile and the target separately. All 
other effects, which would influence the ycM fluctuations 
are neglected. In this way the model gives a lower limit 
for the fluctuations of the initial state CM rapidity. 

From the point of view of global collective flow phe- 
nomena, we would have to consider an initial system 
of particles in local thermal equilibrium. This system 
does not contain non-thermalized, pre-equilibrium emit- 
ted particles, jets, high energy direct gammas, etc. In 
the present estimate we neglect all these effects, as the 
quantitative theoretical estimate of all these effects is ex- 
ceedingly difficult, and even the definition of which parti- 
cles could be considered belonging to the collective initial 
state is not settled. These channels take away consider- 
able energy and momentum from the collective initial 
state, so the center of mass rapidity of the collective ini- 
tial state will be bigger than the "lower limit" estimate 
provided by the model PACIAE. 



A. Particle number asymmetries in PACIAE model 

We first estimate the probability distribution of the 
participant nucleons suffered at least one nucleon-nucleon 
collision. Let us have Na participant nucleons from the 
projectile and Nb from the target. When Na = Nb the 
participant matter is symmetric, so the CM momentum 
and the CM-rapidity vanish. 

At a given impact parameter we have a possibility for 
symmetric fluctuations when Na — Nb change by equal 
number of nucleons. This will not effect the Center of 
Mass. If we have an asymmetry, SN = Na — Nb, this 
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FIG. 2: Initial state fluctuation of the number of extra nu- 
cleons, SN, in 1.38+1.38 A-TeV 0-5% central and 60-70% pe- 
ripheral Pb+Pb collisions in PACIAE model. 



leads to a change of the CM-rapidity. 

Taking into account the effect of overlap geometry and 
of the nucleon-nucleon cross section, the PACIAE model 
fl3| . estimates the SN distribution from N par t fluctua- 
tions as presented in Figures [1] and [5] 

In our model calculations the centrality bins are de- 
fined in terms of the geometrical cross section, b^ ax ir = 
(2Ra) 2 tt, and for example a centrality bin of 60-70% cor- 
responds to an impact parameter range bj], such that 
(bh)/(b 2 max 7r) = 0.6 and (&f 7r)/(6^7r) = 0.7. 

As shown in Figure 1 for the central RHIC collisions 
\SN\/ < Np ar t >~ 1.5%, while for peripheral collisions 
it is 5 %. In peripheral collisions the longitudinally mov- 
ing uppermost and lowermost layers have relatively more 
particles than in central collisions, and so the random 
fluctuations have include relatively more particles, al- 
though the absolute number of particle asymmetry is less. 

In Figure 2 the same results for LHC collisions are 
1.4% for central and 6 % for peripheral collisions. Thus, 
the relative number fluctuation for central collisions de- 
creased slightly due to the difference in the number of 
participants, while for peripheral collisions the small in- 
crease is primarily caused by the difference in the central- 
ity bin. The small difference indicates that the relative 
number fluctuation in peripheral collisions is less sensi- 
tive to centrality bin selection than the absolute numbers. 

At higher energy the cross sections are bigger, so both 
the number of realized primary-primary collisions and 
primary-secondary collisions are bigger. This results in 
an increase in the participant number in the same over- 
lap domain. This leads to the observed fact that while 
the absolute numbers are increasing the relative number 
fluctuations show a smaller increase. 



total momentum of the symmetric part, (Na + Nb— \SN\), 
of the participant matter vanishes. We assume a fix im- 
pact parameter, b and neglect mass number fluctuations 
of the symmetric part of participant matter. Then we 
can assume the mass number of the symmetric part to 
be < Np ar t > — < \SN\ >. If we assume further that 
all of the reaction energy is absorbed in the participant 
matter and < N par t >^> SN then we get 



^VCm(SN) Kt artanh 



SN 



< Npart > 



tanh(yo) 



Thus, the CM-rapidity distribution becomes a series of 
delta functions according to the <5A^-distribution. If we 
allow for the fluctuation of the symmetric mass number 
for a range of impact parameters or a range of multiplici- 
ties, or we allow other channels mentioned above leaking 
energy from the initial state the peaks of the CM-rapidity 
distribution will be smoothed out. 

Figure [3] shows this delta function structure in the re- 
sulting partonic initial state generated by PACIAE model 
for 1.38+1.38 A-TeV 0-5% central Pb+Pb collisions. 
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FIG. 3: Initial state CM-rapidity fluctuation in 1.38+1.38 
A-TeV 0-5% central Pb+Pb collisions in PACIAE model. The 
figure shows that the rapidity change caused by SN—1, 2, 
3, extra nucleons leads to very sharp peak in the CM- 
rapidity distribution, each peak corresponding to a given 5N 
value. This is because there is no tightly bound system to 
absorb energy and momentum in the model. If not so, the 
bound system will allow for rapidity fluctuations at given SN, 
making each sharp peak much wider, and increases the width 
of the overall ycM distribution. 



B. Rapidity fluctuations in PACIAE model 

Let us make a simple estimate: what is the resulting 
CM-rapidity fluctuation. The extra nucleons, 5N, carry 
a longitudinal momentum of Sp z = SN tun sinh(j/o)- The 



The sharp peak structure indicates that all other chan- 
nels (pre-equilibrium emissions, etc.) are neglected in our 
estimate, so the source of rapidity fluctuations is the mo- 
mentum of those extra nucleons, which are not matched 
in originating from the projectile and the target. 
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C. CM-fluctuations of different matter components 

In the partonic initial state generated by the PACIAE 
model a large part of reaction energy is invested into glu- 
ons. The gluons are treated as classical point like parti- 
cles just as the quarks and anti-quarks. If these gluons 
were regarded as a distinct gluon field, then this gluon 
field might keep the partonic initial state system more 
bound and uniform. Then the remaining part (quarks 
and anti-quarks) of the partonic initial state fluctuates 
stronger. 

There are other possibilities, which may increase the 
CM-rapidity fluctuation, e.g. pre-equilibrium emission 
of high energy particles reducing the energy or mass of 
the initial state system; considerable kinetic energy in 
rotation of the initial state system; etc.. 

Gluons have an important role in developing collective 
flow still in the QGP phase (indicated by the constituent 
quark number scaling observed at RHIC). This collective 
flow at high energies may lead to a collective rotation 
[1, [Tcj where a significant part of the collision energy 
remains in longitudinal flow, and so it does not contribute 
to the transverse mass of the system. This would lead to 
a form of collective energy from the gluons, which leads to 
increased ycM fluctuations because this energy reduces 
the transverse mass of the system. Such, collective effects 
are not included in PACIAE, as gluons are treated as 
classical point like particles. 

The initial state fluctuations of the energetic partonic 
matter may be important because the developments of 
these components may not be identical, especially at 
the final FO and hadronization stages of the reaction. 
The gluon fields may contribute to forming the final rest 
masses of the hadrons, and they may contribute different 
amount of thermal and collective kinetic energy to dif- 
ferent hadrons [16] . All effects mentioned above would 
increase the center of mass rapidity fluctuation of the 
initial state, but these are not included in the PACIAE 
model we used. 

Figure @] gives CM-rapidity fluctuation of the quarks 
and anti-quarks in the partonic initial state calculated for 
1.38+1.38 A TeV, 0-5% central and 60-70% peripheral 
Pb+Pb collisions by PACIAE model. The fact that, the 
massive gluon field may carry energy and momentum, 
makes it possible to incorporate part of the fluctuations. 
This enables the model to achieve around a few times 
larger CM-rapidity fluctuations than without a flexibly 
moving massive gluon field as one can see in comparing 
Fig. [3] with Fig. |4j Figure [5] gives the fluctuation of the 
CM-longitudinal momentum per participant nucleon of 
the quarks and anti-quarks in the partonic initial state, 
i.e. p z fluctuation. 

In the PACIAE model calculations above, nearly 57.6% 
of the total collision energy is shared by the quarks and 
anti-quarks and 42.4% by the gluons in the 60-70% cen- 
trality Pb+Pb collisions. These values are 57.9% and 
42.1% for quarks and anti-quarks and gluons, respec- 
tively, in the 0-5% central Pb+Pb collisions. So, how 
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FIG. 4: The CM-rapidity fluctuation of quarks and anti- 
quarks in the initial state calculated for 1.38+1.38 A-TeV, 
0-5% central and 60-70% peripheral Pb+Pb collisions by PA- 
CIAE model. 



gluons are treated is an important issue. 
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FIG. 5: The fluctuation of the CM-longitudinal momentum 
per participant nucleon of the quarks and anti-quarks in par- 
tonic initial state, i.e. p z fluctuation calculated for 1.38+1.38 
A-TeV 0-5% central and 60-70% peripheral Pb+Pb collisions 
by the PACIAE model. 



IV. CONCLUSIONS 

Initial state fluctuations were analyzed in the PACIAE 
model, with particular attention to the CM-rapidity fluc- 
tuations. It was found that in central collisions the longi- 
tudinal asymmetry, arising from different number of pro- 
jectile and target participants, in longitudinal momen- 
tum is around 1.5% only, while for peripheral reactions 
it can reach ±5 — 6% (see Figs. [T]and[2|). In central col- 
lisions the CM-rapidity fluctuations arise from this lon- 
gitudinal asymmetry is not large in the PACIAE model 
as indicated by Figure [31 

We can see in Fig. 0]that the arising CM-rapidity fluc- 
tuation is around ±0.03 units for central collisions and 
around ±0.1 units in peripheral ones. These are about 
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5 — 10 times smaller than the assumptions used in ref. 
j.10] , and this would result in less reduction of the original 
v\(y) calculations. On the other hand, the PACIAE es- 
timates can be considered as conservative lower limits of 
the ycM fluctuations, so the measured ycM fluctuations 
may exceed these values. In the present formulation of 
PACIAE, with point like gluons, the gluon contribution 
would decrease the CM-rapidity fluctuations (cf. Fig. [3]) . 

In the PACIAE partonic initial state study above, we 
do not include the pre-equilibrium emission, the collec- 
tive effects as e.g. rotation, and the formation of excited 
intermediate states. These could lead to the increase of 
CM-rapidity fluctuations. The developing collective flow 
may increase and decrease fluctuations, depending on the 
quantitative details of the developing flow pattern. The 
structure of the collective flow will be detected at the end 
of the reaction, but this pattern develops from the initial 
state in the QGP phase where the gluon component is es- 
sential. The collective flow has both transverse and lon- 
gitudinal components. The pre-collision initial state has 
exclusively longitudinal collective motion. At the time 
point of strongest stopping, this longitudinal flow energy 
is reduced to about 30% of the initial value, while in aver- 
age at the end of the reaction the longitudinal and trans- 
verse energy has about 50 — 50% share 17] . Soft EoS (like 
QGP) and collective rotation may increase the share of 
longitudinal flow energy. The increased longitudinal en- 
ergy (especially from rotation) and the projectile/target 
participant asymmetry may in itself contribute to direct 
increased longitudinal fluctuation. 

The share of longitudinal and transverse flow energies 
also influence the transverse mass of the system, which 
indirectly contributes to longitudinal fluctuations. The 
transverse part of the flow energy increases the trans- 
verse mass, while the longitudinal part reduces it. Larger 
transverse mass reduces the ycM fluctuations. We know 
that with increasing beam energy the collective flow be- 
came more energetic and it is the most dominant phe- 



nomenon at LHC energies. This arises from the initial 
energy and momentum distribution, including the gluon 
components, as these are necessary for the development 
of the large collective flow processes. 

The PACIAE model with point like gluons has less 
ability to incorporate these collective flow effects, and 
about two-thirds of the available energy will contribute 
to the transverse mass, while no direct longitudinal flow 
fluctuation will develop from the initial state asymme- 
tries. Thus, PACIAE with point like gluons underesti- 
mates the ycM fluctuations. 

The initial state longitudinal fluctuations are essential 
for the analysis of the directed flow, as these fluctua- 
tions have significant effect on the measurable wi-flow 
[l0| . The present situation regarding the directed flow is 
rather complex as at RHIC and LHC energies that, the 
observed collective v\ flow is rather weak, \v\\ < 0.001 
at 1] = 0.8, so the vi-flow from the initial state fluctu- 
ations may exceed the global collective v\ flow. Thus, 
the evaluation of vi(pt) at low momenta and low rapidi- 
ties is a complex problem, where the two processes are 
interacting [18[ . The event-by-event longitudinal fluctua- 
tions may be important in the assessment and separation 
of the global directed flow and the directed flow arising 
from the initial state random fluctuations. 



Acknowledgements 

This work was supported by the National Natural Sci- 
ence Foundation of China under Grants No. 11075217, 
11047142, 10975062 and the 111 project of the foreign 
expert bureau of China. L.P. Csernai thanks for the 
kind hospitality of the Institute of Particle Physics of 
the Huazhong Normal University, where part of this work 
was done. 



[1] K. Aamodt et al., (ALICE Collaboration) Phys. Rev. 
Lett. 105, 252302 (2010). 

[2] P.K. Kovtun, D.T. Son and A.O. Starinets, Phys. Rev. 
Lett. 94, 111601 (2005). 

[3] L.P.Csernai, J.I.Kapusta, L.D.McLerran. Phys. Rev. 
Lett. 97, 152303-4 (2006). 

[4] S.A. Bass, A. Dumitru, M. Bleicher, L. Bravina, E. 
Zabrodin, H. Stocker, and W. Greiner, Phys. Rev. C 
60, 021902 (199 9); T. Hirano, P. Huovinen, Y. Nara, 
larXiv:1012.3955V l [nucl-th]; H. Petersen, G.Y. Qin, S.A. 
Bass, and B. Muller, Phys. Rev. C 82, 041901 (2010); B. 
Bauchle and M. Bleicher, Phys. Rev. C81, 044904 (2010). 

[5] H. Petersen, Ch. Coleman-Smith, S.A. Bass, R. Wolpert, 
J. Phys. G 38, 045102 (2011); F.G Gardim, F. Grassi, 
Y. Hama, M. Luzum, J.Y. Ollitrault, Phys. Rev. C 83, 
064901 (2011); P. Bozek, W. Broniowski, J. Moreira, 
Phys. Rev. C83, 034911 (2011); P. Bozek, I. Wyskiel, 
Phys. Rev. C81, 054902 (2010). 



[6] G.Y. Qin, H. Petersen, S.A. Bass and B. Muller, Phys. 

Rev. C 82, 064903 (2010); 
[7] H. Petersen, V. Bhattacharya, S.A. Bass and C. Greiner, 

arXive: nucl-th/1105.0340vl] 
[8] F.G. G ardim, F. Grassi Y. Hama, M. Luzum, J.Y. Olli- 
trault, larXiv:1103.4605l [nucl-th] 
[9] L.P. Csernai, and J.I. Kapusta, Phys. Rev. Lett. 69, 737 
(1992); and L.P. Csernai, and J.I. Kapusta, Phys. Rev. D 
46, 1379 (1992). 
[10] L.P. Csernai, V.K. Magas, H. Stocker, D.D. Strottman, 

arXiv: 1101.3451v3 [nucl-th]. 
[11] V.K. Magas, L.P. Csernai, and D.D. Strottman, Phys. 
Rev. C 64 (2001) 014901; and V.K. Magas, L.P. Csernai, 
and D. Strottman, Nucl. Phys. A 712 (2002) 167. 
[12] L.P. Csernai, Y. Cheng, V.K. Magas, I.N. Mishustin, and 

D. Strottman, Nucl. Phys. A 834, 261c (2010). 
[13] Ben-Hao Sa, Xiao-Mei Li, Shou-Yang Hu, Shou-Ping Li, 
Jing Feng, and Dai-Mei Zhou, Phys. Rev. C 75 054912 



7 



(2007); and Ben-Hao Sa, Dai-Mei Zhou, Yu-Liang Yan, 

Xiao-Mei Li, Seng-Qin Feng, Bao-Guo Dong, and Xu Cai, 

arXiv: 1104.1238vl/nucl-th 
[14] Ben-Hao Sa, A. Bonasera, An Tai, and Dai-Mei Zhou, 

Phys. Lett. B 537, 268 (2002). 
[15] C.B. Chiu, R.C. Hwa, C.B. Yang, Phys. Rev. C 78, 

044903 (2008). 

[16] S. Zschocke, S. Horvat, I.N. Mishustin, and L.P. Csernai, 



Phys. Rev. C 83, 044903 (2011). 
[17] L.Bravina, L.P Csernai, P.Levai, D. Strottman, Phys. 

Rev. C 50, 2161 (1994) 
[18] Contributions of R. Snellings, I. Selyuzhenkov and G. 

Eyyubova et. al. (all for the ALICE collaboration) at the 

Quark Matter 2011 Conf. Annecy France, May 22-28, 

2011, to be published in the proceedings. 



